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Abstract 

o' 

Q^^ ■ Recent results from the CLEO Collaboration on both inclusive and exclusive ra- 

i-^ I diative B decays are used to constrain the parameter space of two versions of the 

Q^' Left-Right Symmetric Model. In the first scenario, when the left- and right-handed 

i-^ ' Cabibbo-Kobayashi-Maskawa mixing matrices are equal, Vl = Vr, the radiative B 

J^ ' decay data is shown to lead to strong bounds on the Wl — Wr mixing angle that are 

rS 
j^ . quite insensitive to either the top quark or Wr mass. The second scenario examined is 

that of Gronau and Wakaizumi wherein 6-quark decays proceed only via right-handed 

currents and Vl and Vr are quite distinct. For this model, the combined constraints 

from Tevatron Wr searches, the B lifetime, and radiative B decays lead to a very 

highly restricted allowed range for the Wl — Wr mixing angle. 
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While the Standard Model(SM) of strong and electroweak interactions is in very good 
agreement with all existing experimental data[|l|, there are many reasons to believe that new 
physics(NP) must exist not far above the scale currently being probed at the SLC, LEP, 
and Tevatron colliders. Although we do not know what form this NP might take there 
are a vast number of proposals in the literature. The best that we can do in the 'pre- 
discovery' era is to use existing data to restrict the properties of this NP and to continue 
searching. While colliders provide us with the capability to directly produce signatures of NP, 
a complementary approach is to hunt for NP indirectly through high precision measurements 
and the observation of rare processes. An excellent working example of such a process has 
been provided us by the CLEO Collaboration which has recently observed the exclusive 
decay B -^ K*'y with a branching fraction of (4.5 ± 1.5 ± 0.9) x 10^^ and has placed an 
upper limit on the inclusive quark-level process of 5(6 -^ 57) < 5.4 x 10^^ at 95% CL. Using 
a conservative estimate for the ratio of exclusive to inclusive decay rates 0, the observation 
of the exclusive process implies the lower bound B{b -^ 57) > 0.60 x 10^^ at the 95% 
CL. These values are, of course, consistent with SM expectations 0] but can used to restrict 
various forms of NP, as has been done in the recent literature]^. It is important to note that 
both the upper as well as the lower bounds can be used to constrain NP since any model 
leading to an extremely suppressed rate for this process is already excluded by the CLEO 
data. Not all of the analyses [^] have taken advantage of this additional constraint. 

One scenario of NP which has been popular in the literature for many years and 
has had many manifestations is the Left-Right Symmetric Model(LRM) p based on the 
extended electroweak gauge group SU{2)l x SU{2)n x f/(l). Amongst other things, this 
model predicts the existence of a heavy, right-handed, charged gauge boson, W^, which can 
in principle mix through an angle with the more conventional W^ present in the SM to 
form the mass eigenstates Wi^2- Data from, e.g., polarized /i decay [0] (in the case of light 



right-handed neutrinos) and universahty requirements H] tell us that the size of this mixing 
must be reasonably small (less than, say, |0| = 0.05, or so) but whose exact magnitude 
depends on the detailed assumptions we make about the other features of the model0. As 
we will see below, the exchange of Wj^ within a penguin diagram, in analogy with the SM W 
exchange, can lead to significant deviations from SM predictions for the b ^ s'y branching 
fraction which is quite sensitive to the both the sign and magnitude of 0. 

In order to numerically determine the branching fraction for b —>■ s'j within the LRM 
there are several sets of parameters whose values we need to address: (i) the mass of the Wr 
itself; (a) the ratio of the right-handed to left-handed SU{2) gauge group coupling constants, 
i.e., K, = Qr/ gi'i {Hi) the mixing angle 0; and {iv) the numerical values for the elements of the 
mixing matrix Vr. For purposes of our discussion below we will treat as a free parameter 
and use the data on the b —>■ s'-f decay itself to constrain as a function of the other degrees 
of freedom, //we assume Vl = Vr, then there are several strong constraints on the Wr mass 
arising from both collider searches] 10] as well as the K^ — Ks mass difference |P, 11] and 



it is likely that Mw^ > l.Qn TeV. Although this possibility is both simple and attractive, 
realistic and phenomenologically viable models can be constructed wherein Vr and Vl are 
quite unrelated as in the scenario of Gronau and Wakaizumi(GW) ]|12| that we will discuss 
in more detail below. In such models, at least some of the conventional constraints on the 
Wr mass can be evaded. However, all such bounds are also dependent on the value of k and. 



within the context of grand unified theories, we generally find that k < l]|T3l- One might 
naively expect in a more general context that this ratio or couplings differ from unity by no 
more than a factor of two or so. 

The approach we follow in performing our calculations has already been discussed in 
our earlier work]!!] and we will refer the interested reader to those papers for calculational 



details. An outline of this approach is as follows. To obtain the 6 ^ 57 branching fraction. 



the inclusive b ^ s'-f rate is scaled to that of the semileptonic decay b -^ Xiu. This 
removes major uncertainties in the calculation associated with (i) an overall factor of ml 
which appears in both expressions and {ii) the various right- and left-handed Cabibbo- 
Kobayashi-Maskawa (CKM) factors. We then make use of the data on the semileptonic 
branching fraction[ll5|, which is given by B{b -^ Xiu) = 0.108, to rescale our result. The 
semileptonic rate is calculated including both phase space(due to the large value of mc/mf,) 



and QCD corrections] 16] with rrih = 5 GeV and rric = 1.5 GeV. The calculation of T{b 



s'j) employs the next-to-leading log evolution equations for the coefficients of the b ^ s 



transition operators in the effective Hamiltonian due to Misiak[17|, the gluon bremsstrahlung 



corrections of Ali and Greub|T8|, the leading corrections from heavy quark effective theory 
(HQET)[jT9[, a running oqed evaluated at the b-quark mass scale, and 3-loop evolution 
of the running a^ matched to the value obtained at the Z scale via a global analysis fll] 
of all data. As we will see, the bounds we obtain on the parameters of the LRM are not 
very sensitive to the remaining uncertainties ||20|| in the calculation of the 6 — > 57 branching 
fraction arising from higher order QCD corrections. In what follows we limit our attention to 
the contributions of the charged gauge bosons to the 6 ^ 57 decay rate. In principle, there 
are potentially other significant contributions in the LRM owing to the extended nature of 
the symmetry breaking sector, i.e., there can be significant contribution from charged Higgs 
exchange as well as from flavor- changing neutral Higgs exchange; we will ignore both these 
possibilities in the analysis below. 

To complete the calculation we use the one-loop matching conditions for the various 
operators [l^ in a form that includes contributions from both the SM and new LRM op- 
erators, i.e., for every 'left-handed' operator present in the SM, the existence of light-right 
symmetry dictates the existence of the corresponding 'right-handed' one. The two sets of 
operators do not mix under QCD evolution and can thus be treated independently. The 



b -^ S'j branching fraction can then be expressed as, using (yQEDirnb) = 132.7, 



^(1 + Q) ^ ^ m + Rj)[{Ll + Rl)f + 2L,i?,^] 



B(h -^ S7) = 7""^^r g(& ^ cfz/)— ^ ' J^ ' ' '7; ' ^ „ F, (1) 



where Q{f,g) is the QCD(phase space) correction to the semileptonic decay b — > ciu. While 
Q as a function of mc/mi, is given in[ffH], the exphcit forms for / and g are given, e.g., in 



/ = {l-y^){l-8y' + y')-24yHny, (2) 

g = -2y[{l-y^){l + 10y^ + y^) + 12y\l + y^)\ny], 

where y = mc/mj,. For y = 0.3, we obtain / ~ 0.520, g ~ —0.236, and Q ^ 2.50(^as(mb)). 
The factor F denotes the relatively small corrections from HQET and gluon bremsstrahlung 
mentioned above and are both of order a few percent. Defining t^ = tancf) and r = 
[MW1IMW2Y (using Mwx — 80.21 GeV in numerical calculations), we obtain for a general 
LRM 

+ \V^n^Hl{l-rf + n\r + tlf], (3) 

2URu{L] + Rl) = 2Kmi~r)'Re{V^,Vi)[{l+rtl) + K\r + tl)]. 

(Note that we implicitly assume that the mass of the right-handed neutrino is sufficiently 
low as to allow its participation in the B decay process and lets us neglect corrections of 
order ml^/m^.) Cj{Jp^ are defined via the low-energy effective Hamiltonian 

-Heff = -^^'^A^rl'PR + C'7{IPL)bF,. , (4) 



where Pl,r = (1 i 75)/2, and whose numerical values are obtained from the operators 
evaluated at the weak scale(~ M^i/J via a renormalization group analysis. This analysis 
is, of course, quite similar to that performed for in the SM case except for the additional 
operators that are present and have non-zero co-efficients at the weak scale. Of course, in 
either the SM or LRM, only a few of these weak scale operator co-efficients are non-zero 
to one loop order. Assuming that the top(t)-quark contribution dominates the penguin 
diagrams (as will be the case in the scenarios we examine below), we obtain in the usual 
notation 



C2r{Mwi 

Ciol{Mwi 
Cior{Mwi 



(l + rtJ)(K6V;:)L, 



R: 



m, 



rrib 
C,ol{Mw,){L ^ R) , 



(5) 



—f^UiV,^,V,f)[A2ix,) - rA^ix^)] + K'iVt,V*)n[tlA,ix,) + rA.ix^)] 



where a;i^2 = ''TT-t/^Wi2- '^^^ co-efficients of the operators corresponding to the gluon pen- 
guin, Csl,r{M'iv-^), can be expressed in a manner similar to Cjl^r^Mw-^^) but with Ai -^ Bf, 
note that both Ai and Bi are the same functions found in the usual SM calculation. Explic- 
itly, we find 
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1 



{x-iy 



^^ I d '^ 2 '^ 2i 

Qt\ -X X H — X H 1 — X ma; 

^H4 2 4 22 



+[l-'+r' 
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3 3, ' 

-X m X 
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(6) 



M^) = li^-^y 



Qtl-^x''-^x + 2 + 3xlnx\ (7) 



'-X + 6x X + 2 — 3x liix 



where Qt = 2/3 is the top-quark electric charge and 51,2(2;) are given by the terms propor- 
tional to Qt in ^1,2(3^)- An important feature to note in the expressions above is the chiral 
enhancement, by a factor of rrit/mb ~ 30, of the terms which involve mixing between the 
Wl and Wr gauge bosons which are proportional to a factor of t,^ = tancj). This implies that 
the decay rate for b ^ s'j should be quite sensitive to small values of t^ even when the W2 
is quite massive. 

We note in passing that the assumption of top-quark dominance of the penguin dia- 
grams may not always be valid in a general LRM since, in principal, the values of the elements 
of both Vl and Vr may conspire to suppress this contribution. This happens, however, in 
only a very small region of the parameter space since rrit/mc > 100. 

Let us first consider the situation where Vl = Vr, in this case the implied lower 
bound on the b ^ s'y branching fraction plays no role in restricting the LRM parameters. 
If we assume that k, = 1 and Mw^ is large, we can ask for the bound on t^ as a function 
of rrit that results from the CLEO limits; this is shown in Fig. 1 for a Wr of mass 1.6 TeV 
and which explicitly displays the b —>■ s'-f branching fraction as a function of t^. Here we 
see that (i) the constraint on the value of t^ is relatively insensitive to rrit and, at 95% CL, 
lies in the approximate range —0.02 < t^ < 0.005. These bounds are much more restrictive 
than what one obtains from either /x decay data (—0.056 < t^ < 0.040) [0 or universality 
arguments (—0.065 < t^ < 0.065) p. (ii) For top masses larger than 120 GeV the b —<■ S'j 
branching fraction(i?) is always found to be in excess of 1.4 x 10~^. These results are found 
to be quite insensitive to the particular values chosen for either the Wr mass or k so long as 



the Wn is reasonably heavy. One may wonder if in fact we can turn this argument around 
in order to get a constraint on M^y^ itself from the CLEO data. To address this issue, we 
fix TUt = 160 GeV with k, = 1 and display B for various values of Mw^ as a function of t^ 
as shown in Fig. 2. Here we see that B itself is not very sensitive to the Wr mass for fixed 
rrit so that no limit is obtainable from this decay mode. Lastly, for fixed rrit = 160 GeV and 
Mwn = 1.6 TeV, we can explore the sensitivity of the resulting bounds on t^ as k, is varied; 
this is shown in Fig. 3, for 0.6 < k < 2. As might be expected the bound strengthens with 
increasing values of k, but only weakly so for positive values of t^. The strengthening of the 
bounds for negative t^ is much more noticeable. It is clear from these figures that the CLEO 
results provide an additional important constraint on the LRM parameters when Vl = Vr is 
assumed and that QCD uncertainties at the level of 10 — 20% will not significantly influence 
the results we have obtained. 



Let us now turn to the perhaps more interesting scenario of Gronau and Wakaizumi 
(GW) wherein B decays proceed only via the right-handed currents. For concreteness we 



take the forms of Vl and Vr as they appear in the original work of GW[12|: 



Vr 



^ 1 A 0^ 

-A 1 

1 



V, 



R 



C —cs 

s{l-c) c^+s2 



V2 



V2 V2 
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\ V2 V2 V2 J 



where A(~ 0.22) is the Cabibbo angle and s ~ 0.09. In order to satisfy B lifetime constraints. 



the parameters in the GW model must satisfy the additional requirement 



Mwi, < 416.2 K 



M 

v^ 



ni/2 

GeV ~ 415 K GeV , (9) 



which arises from recent determinations of Vcb in the SM[^. In addition, to satisfy /x decay 
data, the right-handed neutrino must be sufficiently massive (~ 17 MeV) but this has little 
effect on the B decay itself. Of course, a Wr satisfying the above constraint is relatively light 
and should have a significant production cross section at the Tevatron given the form of Vr. 
In our earlier work we showed that a Wr in the GW model can satisfy both the low energy 
and collider bounds provided that k > 1.5 and Mr > 600 GeV|10| if we assume that the 



Wr decays only into the known SM particles as well as the right-handed neutrino. We will 
respect these conditions when considering the predictions of this model for the 6 — > S7 decay, 
but we should remember that these collider-based limits are clearly softened if additional 
decay modes of the Wr are allowed. 

First, let us fix both M^^/^ and k, in order to satisfy the above constraints and examine 
the predicted value of B in the GW model as a function of t^; this is shown in Figs. 4 and 5 for 
various values of m^. From these two figures we learn that {i) the allowed range of t^f, is more 
tightly restricted in comparison to the Vl = Vr case by both the upper and lower CLEO limits. 
[a) The bounds are quite insensitive to rrit and [iii) the value t^ = is almost excluded by 
the CLEO lower limit. Further, we see as M^^ is increased (also increasing n to satisfy the 
constraints above) we see that the curves become steeper and this forces the allowed regions 
of t(j, to become quite pinched and narrow. In the M^^ = 600(800) GeV case the allowed 
ranges for t^ are found to be -0.43 x 10"^ < t^ < and 0.40 x 10~^ < t^ < 0.81 x 10"^( 
-0.32 X 10^3 < t<,i < and 0.29 x 10"^ < t<,i < 0.60 x 10"^). To say the least, these ranges 
are highly restrictive and it is clear that a more precise determination of the value of B may 
rule out the model as it now stands. To show just how pinched these curves become with 

9 



increasing Mwj^, we fix rrit = 160 GeV and let Mwj^ = 400k GeV while varying k; this is 
shown in Fig. 6. Clearly, as M^y^ grows, the allowed ranges become extremely tight and only 
a very fine tuning of the parameters will allow the GW model to remain phenomenologically 
viable unless other sources of new physics are introduced. 

It is, of course, possible that a modified version of the GW scheme may be realized by 
slightly different versions of both Vl and Vr and several such scenarios exist in the literature. 



Hou and Wyler|]23| have, in fact, two distinct versions of these matrices, denoted by I and II. 
The resulting predictions for the b -^ s'-f branching fraction, B, in both scenarios are quite 
similar and are shown in Figs. 7a and 7b for k = 1.5 and Mwj^ = 600 GeV. (The collider 
bounds on the Wr in both these scenarios are essentially identical to the original GW model.) 
Qualitatively, these predictions are very similar to those of the original GW scheme. Quite 



recently, Hattori et al. have proposed another possible version of the GW scenariol^^ leading 
to the predictions for B in Fig. 8, again assuming k = 1.5 and M\y^ = 600 GeV. In this 
model the 'no-mixing' possibility, t<^ = 0, is completely excluded by the CLEO data, but 
otherwise the results are similar to that of the original GW model. It would seem that a 
general result of the GW approach is to restrict t,^ to very small, but most likely non-zero, 
values. As in the standard GW case, the t^ dependence in both the Hou and Wyler as well 
as the Hattori et al. models becomes somewhat stronger as the Wr mass is increased to 800 
GeV and k is set to 2. 

In this paper we have examined the predictions of the Left-Right Symmetric Model 
for the 6 — > S7 branching fraction in the limit where only the Wj^ and Wj^ gauge bosons 
contribute to the penguin amplitudes. We examined two specific versions of this model, the 
first, wherein left-right symmetry is explicit and Vl = Vr, and the second, in which the 
6-quark essentially decays only through right-handed currents. This corresponds to models 
of the kind first constructed by Gronau and Wakaizumi. In the Vl = Vr case, the limits 
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we obtained on the Wl — Wr mixing angle, 0, were found to be relatively independent of 
the top quark mass and the assumed value of My^^ provided k = 1. For fixed top and 
Wr masses, however, the sensitivity of these constraints to variations in k was found to 
be significant. The bounds we obtained on were comparable to, yet somewhat better 
than, those obtainable from fi decay data or universality arguments. No limit on M^^ is 
obtained from these considerations alone and only the CLEO upper bound was needed to 
obtain the resulting constraints. In the GW-type scenarios, both upper and lower bounds 
on the 6 — i> S7 branching fraction provided important input and were folded together with 
additional constraints arising from Tevatron collider searches as well as the B lifetime. Again, 
the resulting limits on (p were relatively mj-independent and extremely tight, falling into two 
distinct regions in all the cases we examined. An improvement in the CLEO bounds could 
conceivably rule out this scenario if our assumptions remain valid, except, perhaps, for some 
extremely fine-tuned cases. Additional penguin contributions in the form of, e.g., Higgs 
bosons, would then be needed to rescue this approach. 

Perhaps rare B decays may yet provide us with a signature for new physics beyond 
the Standard Model. 
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Figure Captions 

Figure 1. b ^ s'-f decay mode in tlie LRM assuming Vl = Vr as a function of the tangent of the 
Wl — Wji mixing angle, t^. Here we assume k, = 1 and a Wji mass of 1.6 TeV for top- 
quark masses of mt=120(dots), 140 (dashes), 160(dash-dots), 180(solid), or 200(square 
dots) GeV. The horizontal solid line is the CLEO upper bound. 

Figure 2. Same as Fig. 1 but with rrit = 160 GeV held fixed and M^^ varied between 300 (lower 
curve) and 1000 (upper curve) GeV. 

Figure 3. Same as Fig. 1 with rrit = 160 GeV and Mvi/^=1.6 TeV with k, varying between 0.6(left- 
most dotted curve) and 2 (inner- most dash dotted curve). 

Figure 4. Predicted values for the branching fraction(i?) of the b ^ s'y decay mode in the 
Gronau-Wakaizumi version of the LRM as a function of the tangent of the Wl — Wr 
mixing angle, t^. In this figure, k = 1.5 and Mwj^ = 600 GeV is assumed and the 
outer(inner)-most solid line corresponds to rrit = 120(200) GeV and is increased in 
each case by steps of 20 GeV. The dashed horizontal lines are the CLEO upper and 
lower bounds. 

Figure 5. Same as Fig. 4 but with k, = 2 and M^n = 800 GeV. 

Figure 6. Same as Fig. 4 but for rrit fixed at 160 GeV and M^^ = 400/t GeV. Here k is varied 
between 1 and 2 in steps of 0.2 with k = 1(2) corresponding to the outer(inner) curve. 

Figure 7. Predicted values for B in the Hou-Wyler parameterization of Vl and Vr assuming 
K = 1.5 and M^i/^ = 600 GeV. The dotted(dashed, dash-dotted, solid, square-dotted) 
curve corresponds to rrit = 120(140, 160, 180,200) GeV: (a)version I and (b)version II. 

Figure 8. Same as Fig. 7, but for the parameterization of Hattori et al.. 
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